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ABSTRACT 


Initial  steps  are  reported  in  the  development  of  a method  for  deter- 
mining the  instantaneous  velocity  vector  at  the  exit  of  a compressor  rotor 
using  two  Kulite  pressure  probes  with  synchronized  sampling.  A method 
to  determine  the  yaw  angle  of  the  flow  at  any  point  in  the  reference  frame 
of  the  rotor  using  a single  cylindrical  Kulite  impact  probe  and  an  equivalent 
pneumatic  probe  is  described.  The  method  was  verified  in  the  steady  flow 
of  a free  jet  and  applied  in  the  unsteady  flow  field  of  the  compressor.  On- 
line calibration  of  the  Kulite  probe  with  reference  to  the  pneumatic  probe 
was  essential  to  the  success  of  the  method.  Results  for  the  distribution 
of  the  yaw  angle  peripherally  across  the  exit  of  chosen  pairs  of  rotor  blade 
passages,  are  given.  The  analytical  approximation  of  impact  probe 
characteristics,  on  which  the  method  is  based,  and  its  application  with 
experimental  data  are  described. 
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I.  INTRODUCTION 


The  work  reported  here  is  part  of  an  on-going  effort  at  the  Naval 
Postgraduate  School  to  determine  the  aerodynamic  characteristics  and 
performance  of  a single  stage,  axial  transonic  compressor  (Fig.  1)  using 
real  time  instrumentation. 

In  work  already  reported  a multiple  sensor  pneumatic  probe  (Dodge 
probe)  was  developed  to  determine  the  average  velocity  at  the  rotor  exit 
[Ref.  l].  The  design  of  the  probe  was  based  on  a knowledge  of  the  charac- 
teristics  of  simple  impact  probes,  and  experience  which  had  been  gained 
in  representing  analytically  the  characteristics  of  multiple  sensor  probes 
[Refs.  2 & 3],  A new  method  of  representing  the  characteristics  of  mul- 
tiple sensor  probes  was  developed  and  the  time  averaged  flow  from  the 
rotor  was  measured  [Ref.  4],  In  order  to  obtain  information  on  the  flow 
within  the  rotor  itself,  real-time  measurements  were  required.  A syn- 
chronized sampling  system  was  therefore  developed  which  allowed  program- 
mable digital  data  acquisition  from  fixed  instrumentation.  [Ref.  5]  Using 
this  electronic  "pacer,  " sampling  and  analog  to  digital  conversion  of  trans- 
ducer signals  was  controlled  to  be  at  any  of  128  positions  between  any 
selected  adjacent  pair  of  the  18  blades  of  the  rotor,  independent  of  the 
speed. 

A method  of  velocity  measurement  was  proposed  which  used  two  Kulite 
probes  with  synchronized  sampling.  The  technique  would  establish  the 
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instantaneous  and  some  time  dependent  properties  of  the  velocity  field  in 
the  frame  of  the  compressor  rotor.  [Ref.  6]  In  addition  it  was  proposed 
to  relate  the  results  obtained  with  this  method  to  measurements  made 
using  hot  wire  probes  and  using  laser  velocimetry.  The  primary  objec- 
tive of  the  overall  research  effort  was  to  understand  and  to  interpret  the 
phenomena  occurring  at  different  machine  conditions.  [Ref.  6] 

The  purpose  of  the  present  work  was  to  begin  the  development  of  the 
proposed  probe  system  to  measure  "instantaneous"  velocity  at  the  exit 
of  the  transonic  compressor  rotor.  The  ideas  behind  the  two -probe  tech- 
nique were  derived  from  the  previous  work  on  multiple  sensor  probes 
[Ref.  1,2,  3 &4]  and  the  application  of  synchronized  sampling.  Concept- 
ually the  two  probes  can  be  considered  to  supply  measurements  which 
correspond  to  those  of  the  sensors  of  a multiple  sensor  probe.  The 
synchronized  sampling  allows  measurements  from  two  physically  separated 
probes  to  be  taken  (at  different  times)  at  the  same  point  in  the  reference 
frame  of  the  rotor. 

The  initial  step  in  the  proposed  method  was  to  determine  the  yaw 
angle  of  the  flow  at  any  required  point  in  the  rotor  frame.  In  the  present 
work  a method  for  determining  the  flow  yaw  angle  at  a point  in  the  rotor 
frame  using  a single  Kulite  probe  was  developed  and  successfully  tested 
in  the  transonic  compressor. 

Section  II  describes  the  methods  used  to  obtain  velocity  from  pressure 
measurements  taken  with  a multiple  sensor  probe,  and  introduces  the 


10 


proposed  velocity  measurement  method  and  the  method  to  determine  flow 
yaw  angle.  In  Section  III  the  Kulite  and  pneumatic  probes  and  the  data 
acquisition  system  used  in  the  current  work  are  described.  A technique 
for  calibrating  the  Kulite  probe  on-line  to  an  equivalent  pneumatic  probe 
was  developed  and  the  method  for  yaw  angle  measurement  was  verified  in 
experiments  using  a steady,  calibrated  free  jet.  The  tests  and  results 
are  presented  in  Section  IV.  Section  V describes  the  transonic  compres- 
sor test  conducted  to  determine  the  flow  yaw  angle  at  selected  points  in 
the  rotor  frame  and  discusses  the  results.  Conclusions  of  the  study  are 
given  in  Section  VI. 

The  use  of  an  impact  probe  to  measure  yaw  angle  required  that  the 
characteristics  of  such  a probe  be  known.  Measurements  of  the  character- 
istics of  an  impact  probe  were  made  in  a free  jet,  and  an  analytical  expres- 
sion was  found  to  represent  the  results.  The  analytical  representation 
was  fundamental  to  the  successful  measurement  of  the  yaw  angle,  which 
followed.  The  measurements  and  analysis  of  the  results  from  tests  of  a 
simple  impact  probe  are  given  in  Appendix  A. 

The  derivation  of  a numerical  procedure  to  approximate  cylindrical 
probe  characteristics  with  the  empirical  expression  is  given  in  Appendix  B. 
Also  described  is  the  method  used  to  derive  yaw  angle  from  probe  data 
using  the  same  empirical  expression.  Finally,  the  empirical  expression 
was  used  to  derive  the  calibration  characteristics  of  the  Dodge  probe 
which  was  constructed  using  tubing  similar  to  that  used  in  the  construction 


II.  APPROACH 


A.  USE  OF  MULTIPLE  SENSOR  PROBES  TO  MEASURE  VELOCITY 
The  United  Sensor  five  hole  probe  and  the  probe  constructed  by 
F.  J.  Dodge  [Ref.  i],  shown  in  Figure  2 can  be  used  to  measure  air 
average  velocity  when  the  outputs  of  the  multiple  sensors  oriented  at 
various  angles  to  the  probe  axis  are  calibrated  over  a range  of  Mach 
numbers  and  pitch  angles.  Yaw  angle  is  first  determined  directly  by 
balancing  the  outputs  of  the  two  sensors  which  are  set  at  equal  and  op- 
posite angles  to  the  plane  of  the  impact  sensor  and  the  probe  shaft.  Pitch 
angle  (0)  and  non-dimensional  velocity  magnitude  (x)  are  determined  from 
readings  taken  after  the  probe  has  been  rotated  to  balance  the  pressures 

at  sensors  P^  and  P^  (Fig.  2),  and  the  yaw  angle  has  been  noted  from  a 
$ 

vernier  scale. 

The  first  method  for  calibrating  and  applying  multiple  sensor 

probes  uses  the  definitions: 

ft  * ft  (x.  0)  = (P{  - P^)  ^ ^1  (United  Sensor  and  Dodge) 

7 = 7 <x*  = <P1  - P4>  1 <P1  - P22>  (Dodge) 

or 

7=  7 <*.  0)  = (P4  - P5)  / (Pj  - P23>  (United  Sensor) 

♦ Footnote:  The  non-dimensional  velocity  x is  defined  as  x = V/Vt,  where 

V = velocity  magnitude  and  Vt  = "limiting  velocity"  = 2CpT 
for  a perfect  gas,  where  Cp  = specific  heat  at  constant 
pressure  and  T(  = stagnation  temperature. 
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In  the  calibration,  data  for  ft  and  7 are  measured  in  a controlled 
flow  over  the  range  of  pitch  angles  and  Mach  numbers  expected  in  the  un- 
known flow  to  be  measured.  The  data  is  then  approximated  to  yield  fifth 
order  polynomial  expressions  for,  7 - 7 (0,  x)  and  x = x(  ft  , 0).  In 
application,  probe  measurements  in  an  unknown  flow  are  reduced  to  the 
velocity  vector  by  solving  the  polynominals  simultaneously.  This  method 
is  described  in  Refs.  1,  2,  and  3. 

A second  calibration  method  [Ref.  4]  has  been  developed  in  which  (2 
and  7 are  defined  as  in  the  first  method.  Calibration  data  is  obtained  in 
a known  flow  and  reduced  to  fifth  order  polynomial  expressions  for: 

0 = 0 ( 7 * x) 

ft  = ft  (0.  *)  « V (x) 
where  \)  (x)  = ( 7 / 7 x < 1 - x2) 

From  these  approximations,  which  include  the  primary  dependence 
on  velocity  explicitly,  readings  of  ft  and  7 in  an  unknown  flow  are  re- 
duced to  velocity  magnitude  and  pitch  angle  as  in  method  2. 

A third  method  uses  an  analytical  expression  for  the  characteristic 
behavior  of  a single  cylindrical  impact  tube  in  yaw  or  pitch  at  different 
Mach  numbers.  The  overall  calibration  for  a multiple  sensor  probe 
constructed  using  similar  tubing  then  can  be  derived  analytically  as  shown 
in  Appendix  D. 
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B.  MEASUREMENT  OF  VELOCITY  USING  2 KULITE  SEMICONDUCTOR 
PROBES  WITH  SYNCHRONIZED  SAMPLING 

Construction  of  sufficiently  small  multiple  sensor  probes  using  Kulite 
semiconductor  sensors  was  not  considered  to  be  practical  due  to  the  small 
probe  size  required  and  the  low  probability  that  a probe  could  be  construc- 
ted without  the  failure  of  at  least  one  sensor.  Moreover,  such  a probe 
could  not  be  used  to  resolve  peripheral  variations  in  the  yaw  angle  down- 
stream of  the  transonic  rotor. 

What  was  proposed  in  Ref.  6 was  to  use  two  separate  impact  probes 
to  perform  the  tasks  of  the  separate  sensors  in  the  multiple  sensor  probe. 

The  two  probes  are  located  at  different  peripheral  but  identical  axial 
and  radial  stations  between  the  rotor  and  stator  blade  rows.  The  two  pro- 
bes are  shown  in  Fig.  3.  The  probes  are  sampled  under  pacer  control 
such  that,  at  the  time  the  two  readings  are  digitized,  the  probes  are  at 
the  same  point  in  the  reference  frame  of  the  rotor.  [Ref.  6]  A computer 
peripheral  device  which  is  required  to  control  the  data  acquisition  (pacer) 

was  developed  by  West.  [Ref.  5] 
o 

The  90  or  type  A probe  shown  in  Figure  3 can  be  rotated  in  yaw 

about  the  tip  to  angles  corresponding  to  the  sensors  for  Pj,  P£,  and  P^ 

in  the  Dodge  probe.  The  55°  or  type  B probe  shown  in  Figure  3 must  be 

set  at  the  balanced  or  zero  yaw  angle  (determined  using  the  type  A probe) 

to  measure  the  pressure  corresponding  to  P in  the  Dodge  probe.  To 

4 

determine  velocity,  one  of  the  three  methods  described  in  the  preceding 
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section  must  be  used  to  represent  and  apply  the  calibration.  The  two 


probes  must  be  calibrated  together  in  a steady  uniform  flow  as  if  they 
were  the  sensors  of  a multiple  sensor  probe. 

The  initial  step  in  applying  the  two  probe  system  to  measure  velocity 
is  to  determine  the  yaw  angle  of  the  flow  at  the  required  point  in  the  rotor 
frame.  Only  then  can  the  two  probes  be  set  to  angles  at  which  the  cali- 
bration was  established.  In  the  present  work,  a method  for  determining  the 
yaw  angle  at  each  point  in  the  rotor  frame  was  developed  using  only  the 
type  A probe. 

Since  the  method  makes  use  of  the  characteristics  of  cylindrical  impact 
probes,  a preliminary  investigation  was  made  of  the  characteristics  of 
pneumatic  cylindrical  impact  probes.  The  results  are  given  in  Appendix  A. 
It  was  found  that  the  response  of  the  probe  could  be  represented  as  a pres- 


sure coefficient  Cp,  defined  as: 

F^p-F*.  .01 


Cp  = 


Ip  M 


AM 


[ sin2  Bit-  ^0)]” 


where 


A,  B and  N are  constants  for  the  particular  probe, 
M = Mach  number 
p = static  pressure 
~jf  - ratio  of  specific  heats 
Pp  = probe  pressure 
= impact  pressure 


16 


V = angle  of  the  flow  to  a reference  scale 

angle  of  the  axis  of  the  probe  to  the  reference  scale 
When  the  flow  angle  to  the  probe  axis  was  a combination  of  pitch  angle, 
0,  and  yaw  angle,  0(  , then  V was  calculated  using  the  geometrical 
relationship: 

y = cos  "*(cos  0 cos (X  ) 

These  results  provided  the  basis  for  the  development  of  a technique 
to  derive  yaw  angle  from  Kulite  probe  measurements:  the  development  is 
described  in  Section  IV,  following  a description  of  the  probes  and  data 
acquisition  method  given  in  the  next  section. 
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m.  INSTRUMENTATION 

A.  PROBE  DESCRIPTION 

Figure  3 is  drawing  of  the  XB-062-25  semiconductor  probes  manufac- 
tured by  Kulite  for  the  tests  conducted  in  this  study.  The  sensors  in  each 
probe  tip  were  standard  CQ-052  series  Kulite  ultra  miniature  pressure 
sensors.  Kulite  Type  "B"  screens  were  installed  on  the  probes  to  protect 
the  sensors  from  particle  impingement  in  high  velocity  flows.  Figure  4 
shows  the  probe  tip  dimensions  and  internal  design. 

Kulite -equivalent  pneumatic  probes  were  constructed  at  Naval  Post- 
graduate School  that  were  similar  (to  within  0.  002  inches)  to  the  external 
geometry  of  the  Kulite  probes  shown  in  Figures  3 and  4.  Initial  impact 
pressure  measurements  taken  with  the  Type  A Kulite  and  pneumatic 
equivalent  probes  were  consistently  10%  below  actual  total  pressure  in 
the  flow.  It  was  immediately  recognized  that  flow  stagnation  was  occurring 
in  the  center  of  the  probe  face  and  not  at  the  holes  in  the  screen.  There- 
fore, a shroud  of  0.083  inch  outside  diameter,  stainless  steel  tubing  was 
added  as  shown  in  Figure  4.  Subsequent  tests  indicated  good  agreement  of 
impact  pressure  measurements  with  actual  total  pressure  with  no  notice- 
able loss  in  the  frequency  response  of  the  Kulite  probe. 

B.  DATA  ACQUISITION  SYSTEM 

Data  acquisition  from  the  Kulite  probes  and  conventional  strain  gauge 
transducers  was  controlled  by  a Hewlet-Packard  21  MX  computer.  A pacer 


I 
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developed  by  West  [Ref.  5]  was  used  for  "synchronized"  sampling  behind 
the  transonic  compressor  rotor  as  shown  in  Figure  5.  Synchronized 
sampling  allows  the  properties  of  the  flow  to  be  determined  at  a point 
which  is  fixed  with  respect  to  the  rotor.  The  pacer  uses  conditioned 
1 per  blade  and  1 per  rev  signals  to  generate  a trigger  pulse  at  a selected 
location  in  the  rotor  frame.  There  are  256  equally  spaced  selectable 
locations  in  each  rotor  blade  pair.  Steady  and  nonsteady  flow  properties 
can  be  determined  by  sampling  at  a discrete  location  as  many  times  as 
is  necessary.  Figure  6 shows  the  complete  data  acquisition  system  used 
in  the  present  study.  Appendix  C is  a discussion  of  the  software  used 
with  the  21  MX  computer. 


IV.  PRELIMINARY  MEASUREMENTS  IN  A FREE  JET 


A.  INTRODUCTION 

In  order  to  use  semiconductor  probes  for  quantitative  pressure  meas- 
urements in  an  environment  in  which  the  temperature  is  unknown,  a means 
of  on  line  calibration  must  be  used.  The  output  voltage  of  the  Kulite  probe 
is  linear  with  respect  to  the  pressure  difference  across  the  sensing 
diaphragm.  A typical  Kulite  calibration  is  shown  in  Figure  8.  However, 
large  shifts  in  output  voltage  due  to  temperature  changes  have  been  re- 
ported and  were  confirmed  by  Paige  [Ref.  7].  For  this  reason,  Kulite 
semiconductors  have  not  been  used  for  absolute  pressure  measurements 
in  turbo  machinery  but  have  been  used  to  measure  the  fluctuations  in  un- 
steady or  periodic  flows.  In  calibration  tests  conducted  in  a steady  air 
flow  early  in  the  present  work,  these  shifts  of  level  of  output  were  ob- 
served. Shifts  in  the  slope  of  the  output  were  also  observed  when  the 
Kulite  probe  was  rotated  to  a new  flow  angle.  Therefore,  a technique  to 
calibrate  the  Kulite  probe  against  a geometrically  identical  pneumatic 
probe,  similarly  oriented  in  the  flow,  using  the  computer  on  line,  was 
developed. 

The  goal  of  the  experiments  described  in  this  section  was  to  develop 
the  on-line  calibration  procedure  and  to  verify,  in  steady  flow,  a method 
for  measuring  yaw  angle  in  the  compressor.  A four  inch  free  jet  was 
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used  to  provide  a known  uniform  flow  field  for  the  tests.  The  test  appara- 
tus  is  described  in  reference  8.  The  Kulite  type  A and  the  equivalent 
pneumatic  type  A probes  were  inserted  in  the  flow  to  be  two  inches  apart 
on  opposite  sides  of  the  jet  axis.  The  mounting  apparatus  allowed  both 
probes  to  be  pitched  and  yawed  simultaneously  and  to  be  set  to  similar 
angles  with  respect  to  the  flow  direction.  A Prandtl  probe  was  used  to 
monitor  the  jet  velocity. 

B.  KULITE- EQUIVALENT  PNEUMATIC  PROBE  TEST 

Tests  in  the  free  jet,  similar  to  those  described  in  Appendix  A,  were 
conducted  to  establish  the  characteristics  of  the  type  A equivalent  pneu- 
matic probe  with  respect  to  yaw  angle  and  Mach  number.  The  range  of 
Mach  numbers  surveyed  was  . 4 to  .6. 

The  results  are  shown  in  Figure  7.  The  pressure  coefficient  vs. 
yaw  angle  characteristics  for  this  probe  and  the  variation  with  Mach 
number  were  found  to  be  qualitatively  similar  to  those  of  the  cylindrical 
impact  probe  reported  in  Appendix  A.  The  coefficients  A,  B,  N,  and 
were  successfully  calculated  as  described  in  Appendix  B.  It  is  concluded 
therefore  that  the  technique  for  calculating  the  zero  yaw  angle  described 
in  Appendix  B could  also  be  used  with  the  type  A probe  geometry. 

C.  KULITE  PROBE  TESTS 

1.  Averaging  Techniques 

The  first  measurements  from  the  Kulite  probe  in  a free  jet  veri- 
fied the  high  level  of  turbulence  and  unsteadiness  which  were  known  to  be 
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present  from  earlier  measurements.  In  order  to  investigate  the  charac- 
teristics of  the  Kulite  probe  with  respect  to  yaw  angle  and  Mach  number, 
some  method  of  averaging  over  a number  of  data  samples  was  necessary. 

The  purpose  of  the  first  test  therefore  was  to  determine  the  effect 
of  sample  number  and  interval  on  the  average  measurement  obtained  from 
the  Kulite  probe. 

At  one  test  condition,  three  ten  minute  tests  were  conducted  in 
which  the  techniques  for  sampling  the  Kulite  probe  were  varied.  During 
the  tests,  the  Kulite  equivalent  pneumatic  probe  pressure  was  read  at 
close  intervals  using  a water  column  U-tube  manometer  referenced  to 
atmosphereic  pressure.  The  Kulite  reference  pressure  was  sensed  by 
a conventional  strain  gauge  transducer  connected  to  and  calibrated  on 
one  input  channel  of  the  A/D  converter.  The  reference  pressure  was 
obtained  at  one  minute  intervals  by  recording  the  average  100  samples 
taken  at  10  millisecond  intervals. 

In  each  of  the  three  tests,  a data  point  for  the  Kulite  probe 
pressure  was  recorded  at  1 minute  intervals  for  a period  of  10  minutes. 
Each  data  point  consisted  of  taking  the  average  voltage  from  a number 
of  ensembles  of  samples  taken  at  10  microsecond  intervals.  The  data 
for  the  methods  used  in  the  three  tests  is  shown  in  the  following  table: 
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SAMPLES  PER 
METHOD  ENSEMBLE 

# OF 

ENSEMBLES 

ENSEMBLE 

INTERVAL 

(SEC) 

TOTAL  TIME 
(SEC.  ) 

1 

100 

1 

- 

1 

2 

1000 

1 

- 

10 

3 

100 

10 

2 

20 

For  each  data  point. 

the  pressure 

was  calculated 

using  the 

average  voltage  from  the  Kulite  probe  and  the  recorded  reference  pres- 
sure  as  shown  in  the  following  section,  C.2.  The  results  are  shown  in 
Figure  9. 

It  can  be  seen  that  methods  two  and  three  gave  acceptable  accuracy 
whereas  method  one  did  not.  The  failure  of  method  one  was  the  result  of 
high  frequency  turbulence  in  the  free  jet.  Seen  on  the  oscilloscope,  the 
dominant  fluctuations  in  the  Kulite  signal  had  a period  of  about  two  milli- 
seconds, so  that  a sample  duration  of  one  millisecond  could  not  give  a 
correct  time  average.  Method  two  was  used  in  subsequent  tests  in  the 
free  jet. 

2.  Calibration  Method 

The  need  to  calibrate  the  Kulite  transducer  to  account  for  changes 
due  to  temperature,  and  for  shifts  in  output  slope  which  occurred  when 
the  probe  was  rotated  in  the  flow,  was  stated  in  Section  IVA.  A calibra- 
tion procedure  was  therefore  developed  in  the  free  jet  tests  and  subsequen- 
tly applied  in  the  measurement  of  yaw  angle.  The  Kulite  probe  calibra- 
tion used  the  assumption  that  the  equivalent  pneumatic  probe  measured 
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the  time  average  of  the  pressure  on  the  face  of  the  Kulite  probe  when  the 
two  probes  were  similarly  oriented  in  the  flow.  The  calibration  procedure 
was  to  apply  a controlled  pressure  to  the  reference  tube  of  the  Kulite 
probe,  then  sample,  (as  described  in  Section  IV  C.  1)  the  pressure  from 
the  equivalent  pneumatic  probe,  Pp,  the  Kulite  reference  pressure,  Pr, 
and  the  Kulite  output  voltage.  This  procedure  was  repeated  for  a number 
of  calibration  points,  changing  the  reference  pressure  for  each  point.  The 
data  was  reduced  on-line  using  a linear  least  squares  routine  to  obtain 
the  constant  coefficients  xQ  and  xj  in  the  equation, 

Pp  - Pr  = x0  + xj  E" 

where  E was  the  average  amplified  output  voltage  of  the  Kulite  transducer. 

In  subsequently  applying  the  calibration,  samples  of  Kulite  output 
voltage  were  reduced  on  line  to  values  of  pressure,  Pk,  using  the  calibra- 
tion coefficients  x0  and  Xj  in  the  equation 
Pk  = xQ  + xj  * E + Pr 

An  experiment  was  conducted  to  determine  the  minimum  number 
of  reference  pressure  settings  required  to  calibrate  the  Kulite  probe  with 
acceptable  accuracy.  A calibration  test  was  carried  out,  as  described 
above,  using  eight  different  reference  pressures.  The  coefficients  xQ  and 
Xj  were  computed  first* using  all  eight  data  points.  Then,  using  the  same 
data,  new  coefficients  were  computed  using  various  subsets  of  the  eight 
data  points.  Table  IV- 1 shows  the  subsets  used  and  the  deviation  of  the 
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results  from  the  eight  point  calibration.  The  results  indicated  that  only 
small  improvements  in  the  accuracy  were  obtained  by  using  more  than 
two  points.  Two  points,  if  properly  chosen  to  include  the  range  of  ex- 
pected pressures,  gave  a sufficiently  accurate  calibration.  In  subsequent 
experiments  a two-point  calibration  was  therefore  used  with  atmospheric 
pressure  as  one  reference  pressure  and  the  flow  total  pressure  as  the 
other. 

D.  YAW  ANGLE  DETERMINATION 

The  purpose  of  this  experiment  was  to  verify,  in  steady  flow,  a method 
of  measuring  the  zero  yaw  angle  behind  the  rotor  of  the  transonic  compres- 
sor using  the  type  A Kulite  probe.  Measurements  were  made  in  the  four 
inch  free  jet  at  Mach  numbers  of  0.406  and  0.  587.  At  each  Mach  number 
both  probes  were  set  at  pitch  angles  of  -10  , 0°,  and  +10  . At  each  pitch 
angle  both  probes  were  set  at  yaw  angles  of  +60  , +45°,  +30  , ^+15  , and 
0°  to  the  flow.  After  each  change  of  angle  the  Kulite  probe  was  calibrated 
as  described  above.  Water  column  readings  of  the  Prandtl  probe  and  the 
Type  A pneumatic  probe  were  recorded,  the  Type  A pneumatic  pressure 
was  sampled  and  the  Type  A Kulite  average  voltage  was  sampled  and 
reduced  on  line  to  values  of  pressure. 

In  the  off-line  analysis,  zero  yaw  angle  was  calculated  from  the  data 
using  the  least  squares  method  described  in  Appendix  B.  The  number  of 
points  needed  to  determine  the  zero  yaw  angle  was  investigated  by  varying 
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the  data  points  included  in  the  least  squares  calculation.  The  results 
are  given  in  Table  IV-2  for  3,  5,  7 and  9 points. 

It  was  seen  that  the  zero  yaw  angle  was  determined  to  within  0.  5 
degrees  for  pitch  angles  of  -10,  0 and  +10  degrees.  It  was  concluded  that 
the  Kulite  probe,  when  carefully  calibrated  to  an  equivalent  pneumatic 
probe,  can  be  used  to  measure  yaw  angles  with  reasonable  accuracy. 

The  need  to  calibrate  the  Kulite  probe  at  each  orientation  was  con- 
firmed by  an  inspection  of  the  calibration  coefficients  determined  in  the 
experiment.  The  coefficients  are  shown  in  Table  IV-3. 

The  variation  in  the  coefficients  appeared  to  be  sufficiently  random 

in  nature  as  to  preclude  any  attempt  to  describe  the  change  as  a function 

\ 

of  Mach  number  and  flow  angle.  The  variations  were  not  understood, 
however,  the  technique  of  recalibrating  at  each  angle  was  one  which  could 
also  be  applied  in  the  compressor.  It  should  be  noted,  however,  that  in 
the  periodic  flow  behind  a compressor  rotor,  application  of  the  same 
calibration  and  measurement  techniques  will  force  agreement  between 
the  absolute  values  of  the  pneumatic  time  average  and  Kulite  time  average 
pressures. 
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(i)  DATA 


POINT 

P-Pr  (ins. 

water) 

AMPLIFIED  KULITE  OUTPUT 

(millivolts) 

l 

-42.  8 

-68.9415 

2 

-27.  5 

-46. 193 

3 

-12.  5 

-25.2698 

4 

- 0.2 

- 6.05109 

5 

24.  2 

29.0726 

6 

36.0 

46. 186 

7 

47.  2 

62.0352 

8 

57.  8 

77.7253 

Impact  pressure  at  the  probe  face. 

pt=37.9  ins.  water  gauge. 

Least  squares  fit  to  linear  calibration  given  by,  P-Pr=xQ+XjE  : 

xq=4.  39565  ins 

. water, 

=686.  358  ins.  water/volt. 

(ii)  Least  Squares  fit  to  linear  calibration,  P-P  =x  +x  E,  using 

subsets  of  points. 

Points  used  for  Subset 

(xo-x0> 

(Ej-x^/xj  (%) 

1.3, 5,7 

0.  2 

0.  16 

1.4,8 

0.  1 

0.  01 

2.4,7 

-0.  1 

-0.  60 

3,4.6 

0.  0 

0.  70 

1.3.5 

0.  1 

0.  44 

3,5.7 

0.  3 

0.  50 

1,8 

0.  2 

0.  07 

2,7 

0.  0 

-0.  51 

3.6 

0.3 

1.  14 

2.5 

-0.  1 

-0.  36 

1.5 

o 

• 

o 

-0.  10 

2,  6 

-0.  1 

0.  39 

3,7 

0.  4 

-0.07 

4.6 

-0.  2 

-0.  11 

Table  IV -1  Results  of  Kulite  probe  calibration  in.  steady  flow 
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P^=Kulite  pressure  calculated  using  linear  calibration 

Pjce=Kulite  equivalent  pneumatic  pressure  measured  with  conventional 
transducer 


(P.  ) =Kulite  equivalent  pneumatic  pressure  measured  with  water 
ke  m , 

column  manometer 

Entries  in  the  table  are  the  zero  yaw  angles  ( (y  Q),  in  degrees,  calculated 
from  the  above  measurements  using  subsets  of  the  data  recorded  at  a total 
of  nine  angles. 

The  subsets  are  as  follows: 

Subset  # Angles  included 

1 0°, +15°, +30°, +45°, +60° 

2 0°, +30°, +45°, +60° 

3 0°, +45°, +60° 

o o 

4 0 , +60 

Subset  No.  of  points  0(  from  P.  (Y  from  P,  (V  from  (P,  ) 

o k o ke  ^o  ke  m 

(M=0.  406,  0=0°) 


1 

9 

.70 

. 80 

.81 

2 

7 

.75 

.80 

.81 

3 

5 

.82 

.86 

. 86 

4 

3 

1.31 

1.29 

1.33 

(M=0.  406,  0=+lO°) 

• 

1 

9 

1.38 

1.  10 

1.08 

2 

7 

1.40 

1.11 

1.09 

3 

5 

1.48 

1.  15 

1.  14 

4 

3 

2.09 

1.61 

1.60 

(M=0.  406,  0 = -lO°) 

1 

9 

.46 

.60 

.60 

2 

7 

.49 

.62 

. 60 

3 

5 

.60 

.65 

.62 

4 

3 

.78 

.98 

.89 

Table  IV-2  Calculated  zero  yaw  angles  from  test  described 

in  Section  IV -D. 
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Subset 

No.  of  points  0£Qfrom 

(M=0.  587,  0=0°) 

CCjTOm  Pke 

0£o£rom  (p 

1 

9 

.86 

.85 

.85 

2 

7 

.86 

.85 

.86 

3 

5 

.95 

.93 

.92 

I 

4 

3 

1.33 

(M=0.587,a  = 10°) 

1.47 

1.44 

9 

.91 

.79 

.83 

2 

7 

.93 

.81 

.84 

3 

5 

.95 

00 

e 

.88 

4 

3 

1.38 

(M=0.  587,  0 = -lO°) 

1.  17 

1.25 

1 

9 

.85 

.77 

. 73 

2 

7 

.85 

.78 

.74 

3 

5 

O' 

00 

• 

.82 

.78 

4 

3 

1.35 

1.28 

1. 23 

IV. 2 (cont)  Calculated  zero  yaw  angles  from  test  described 

in  Section  IV-D. 
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Probe  Yaw 
Angle  (Deg) 

a=o° 

X 

o 

X1 

0=- 

X 

o 

(M=0.  406) 

10° 

X1 

0 = 

X 

o 

+ 10° 

X1 

0 

-3.9 

696 

-3.9 

686 

-5.8 

685 

15 

-5.5 

721 

-4.2 

694 

-5.9 

691 

-15 

-3.7 

686 

-5.2 

691 

-6.3 

694 

30 

-2.6 

683 

-2.8 

693 

-5.  1 

690 

-30 

-5.0 

689 

-5.3 

692 

-7.4 

691 

45 

-0.5 

685 

-0.2 

689 

-2.3 

694 

-45 

-7.0 

686 

-7.4 

698 

-9.4 

677 

60 

+4.5 

715 

+3.  6 

681 

+3.  1 

679 

-60 

-9.3 

711 

-9.6 

693 

-12.4 

671 

(M  = 0.587) 

0 

-4.8 

689 

-5.0 

692 

-4.5 

685 

15 

-4.0 

690 

-5.0 

694 

-5.3 

684 

-15 

-6.  5 

691 

-6.  1 

689 

-6.  1 

687 

30 

-3.5 

692 

-1.3 

697 

-3.8 

681 

-30 

-9.  1 

692 

-8.5 

690 

-10.0 

691 

45 

-2.0 

735 

+5.5 

684 

+0.9 

692 

-45 

-11.  1 

688 

-12.6 

685 

-13.7 

688 

60 

+ 16.  1 

673 

+ 16.7 

684 

+ 16.  3 

672 

-60 

-20.0 

697 

-20.4 

680 

-19.6 

659 

x in  inches  water 
o 

Xj  in  inches  water /volt 


Table  IV-3  Variation  of  Kulite  calibration  coefficients  in 
yaw  angle  test  on  the  free  jet. 
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V.  FLOW  ANGLE  MEASUREMENTS  IN  A TRANSONIC  COMPRESSOR 


r. 


A.  TRANSONIC  COMPRESSOR 

The  transonic  compressor  is  shown  in  Figure  1.  It  is  driven  by  an 
air  turbine  drive  unit  capable  of  supplying  450  horsepower  at  30,000  RPM. 
When  operating  at  the  design  point,  the  relative  Mach  number  at  the  rotor 
blade  tip  is  1.5.  The  flow  rate  is  controlled  by  an  electrohydraulic  rotat- 
ing throttle  plate  located  at  the  inlet  duct,  which  also  contains  a filter  and 
flow  measuring  nozzle.  An  Allis -Chalmers  multistage  axial  compressor 
supplies  the  turbine  drive  air.  A complete  description  of  the  test  facilities 
is  given  in  reference  8. 

B.  PROBES  AND  INSTRUMENTATION 

A Kulite  type  A probe,  a pneumatic  type  A probe  (Fig.  3),  and  the 

Dodge  probe  (Fig.  2)  were  inserted  into  the  compressor  at  an  axial  distance 

of  0.65  inches  (.31  chord  lengths)  behind  the  rotor  trailing  edge.  A similar 

probe  installation  is  shown  in  Figure  10.  The  radial  displacement  of 

each  probe  tip  from  the  case  wall  was  0.  9 inches.  The  three  probes  were 

o 

separated  peripherally  at  45  intervals. 

The  pneumatic  connections  to  the  Dodge  probe  PI  sensor  and  to  the 
pneumatic  type  A probe  are  shown  in  Figure  11.  The  pneumatic  pressures 
were  sensed  by  conventional  differential  strain  gauge  transducers  referenced 
to  atmosphere. 


31 


The  two  strain  gauge  transducers  and  the  Kulite  type  A probe  trans- 
ducer were  connected  through  conditioning  circuits  to  separate  input 
channels  of  the  high  speed  data  system  shown  in  Figure  6.  The  two 
transducers  were  calibrated  using  a water  column  manometer  as  a 
standard.  The  calibration  coefficients  obtained  were  used  in  the  on-line 
program  to  reduce  the  transducer  voltage  outputs  to  pressure  in  inches  of 
water  gauge.  A control  program  was  written  and  used  with  the  program 
"DATACQ"  to  acquire  data  during  the  tests.  These  programs  are  des- 
cribed and  listed  in  Appendix  C.  On-line  calibration  of  the  Kulite  trans- 
ducer was  an  integral  part  of  the  experimental  procedure. 

C.  TEST  PROCEDURE 

The  compressor  was  stabilized  at  the  desired  operating  condition  and 

steady  state  performance  data  were  recorded.  The  Dodge  probe  was 

rotated  to  balance  the  pressures  at  sensors  P2  and  P3  and  the  yaw  angle, 

CC  , was  recorded.  This  angle,  0^  , was  taken  as  a reference  angle 

for  the  Type  A probes.  The  pneumatic  and  Kulite  Type  A probes  were 

rotated  together  to  the  following  angle  settings,  in  turn:  ± 30°  * 

_ o — o 

C<+  45  ,Q(+  55  . At  each  setting,  the  following  procedure  was  carried  out: 

1.  The  Kulite  probe  output  was  sampled  in  Pacer  "Free  Run"  mode. 
1681  samples  were  taken  at  10  microsecond  intervals,  the  aver- 
age was  calculated  and  reduced  on  line  to  pressure. 

2.  The  pressure  at  the  PI  sensor  of  the  Dodge  probe  and  the  pneu- 
matic equivalent  pressures  were  sampled. 
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3.  Kulite  output  was  sampled  in  a survey  at  128  locations  across 
blade  pair  2 using  synchronized  sampling.  Each  location  was 


sampled  10  times,  (on  successive  revolutions)  and  the  average 
of  the  10  samples  was  reduced  on  line  to  pressure. 

4.  The  pressure  survey  data  were  displayed  on  the  HP  9862A  plot- 
ter and  stored  in  the  HP9768  mass  memory  system. 

5.  Steps  3 through  5 were  repeated  for  blade  pair  8. 

Data  recorded  during  the  test  are  given  in  Table  V-l.  Plots  of  the 
data  obtained  by  synchronized  sampling  are  shown  in  Figures  12  and  13. 
The  synchronized  data  was  observed  to  agree  qualitatively  with  oscil- 
loscope traces  of  the  Kulite  output  voltage. 

D.  DATA  REDUCTION  TO  FLOW  ANGLE 

The  data  from  seven  angles  were  used  to  calculate  a zero  yaw  angle 
for  each  discrete  location  within  the  blade  pairs,  using  the  method  des- 
cribed in  Appendix  B.  The  BASIC  program  "KAW78"  used  to  reduce  the 
data  is  listed  in  Table  V-2. 

E.  RESULTS  AND  DISCUSSION 

The  distributions  of  the  zero  yaw  angle  for  the  Kulite  probe  (the  yaw 
angle  of  the  flow  relative  to  the  axial  direction)  for  blade  pairs  2 and  8 
are  shown  in  Figures  14  and  15  respectively. 

In  both  figures,  the  wakes  of  the  rotor  blades  are  evident,  and  fluctua- 
tions are  seen  to  be  present  outside  of  the  wakes.  It  should  be  noted  that 
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the  shape  of  the  observed  angle  variation  through  the  wake  would  be  the 

result,  qualitatively,  of  the  flow  velocity  being  non-uniform  in  magnitude 

but  at  a constant  angle  relative  to  the  rotor.  The  fluctuations  outside 

the  wakes  are  not  understood.  They  are  seen  to  be  more  definite  in  blade 

pair  2 than  in  blade  pair  8.  Further  analysis  of  the  data  and  more  data 

at  different  test  conditions  are  needed  to  determine  the  accuracy  of  the 

detail  which  is  evident  in  Figures  14  and  15. 

The  numerically  averaged  yaw  angles  for  blade  pairs  2 and  8 were 
o o 

22.8  and  22.  5 respectively.  The  pneumatically  averaged  yaw  angle 

o 

obtained  from  the  Dodge  probe  was  23.2  . The  yaw  angle  calculated  using 

o 

time  averaged  pressures  from  the  equivalent  pneumatic  probe  was  22.  2 . 
This  agreement  supports  the  magnitudes  of  the  yaw  angles  shown  in 
Figure  14  and  Figure  15,  for  which  no  other  verification  is  presently 
available. 
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Compgressor  Performance  Data 

N = 18360  RPM 

Weight  Flow  = 11.92  lb/aec 

Ttl  =76*  1<>F 

Ptj  = 389.7  inches  water 
Pressure  Ratio  (t-t)  = 1.  167 
Efficiency  (t-t)  = 0.92 

Data  From  Probe  Measurements 


Pressures  in  inches  of  water 


Point 
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• 

TABLE  V-l  TRANSONIC  COMPRESSOR  TEST  DATA 
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TABLE  V-2.  BASIC  PROGRAM  "KAW781 


200  FOR  1=1  TO  7 

210  LET  C9=C0S (B*(A[ l ,1 ] -YO) ) 

220  LET  C=A12,Il-A*(l-C9t2)tN 

230  LET  81 =2*A*N*(A( 1 f I ] -Y0)*( l-C9T£)T (N  — 1 )*C9*SIN (o*(Al l.IJ-YO)) 
240  LET  T(K] =C*Bl*a/(Al 1 ,1 1 -YO)+TC K 1 


VI.  CONCLUSIONS 


In  measurements  made  downstream  of  the  transonic  compressor 
rotor,  there  was  found  to  be  good  agreement  (better  than  1%)  between  the 
yaw  angle  determined  using  a Dodge  pneumatic  probe  and  the  average  of 
the  peripheral  distribution  of  the  yaw  angle  calculated  from  a system  of 
Kulite  probe  measurements.  It  was  therefore  concluded  that  the  tech- 
nique to  determine  "instantaneous"  yaw  angle  proposed  here  and  verified 
in  a steady  free  jet,  also  gave  good  results  in  the  unsteady  compressor 
flow  field.  The  yaw  angle  variation  in  the  rotor  blade  wakes  was  resolved 
and  the  results  contained  very  little  scatter.  A weak  periodic  structure 
was  detected  across  one  of  the  two  pairs  of  blade  passages  measured 
(Blade  pair  #2).  More  analysis  of  the  data  is  required  however,  before 
conclusions  can  be  drawn  concerning  the  detail  which  appears  to  be  present 
in  the  results. 

The  yaw  angle  measurement  was  a necessary  step  in  the  development 
of  a method  to  determine  the  velocity  field  at  the  rotor  exit  from  stationary 
instrumentation.  The  next  step  is  to  determine  the  best  method  for  repre- 
senting and  applying  the  calibration  of  a two  probe  system.  The  fluctua- 
tions measured  in  the  zero  yaw  angle  (or  turning  angle)  behind  the  trans- 
onic compressor  rotor  makes  less  attractive,  but  does  not  preclude,  a 
calibration  method  that  requires  the  type  B probe  to  be  set  at  the  zero 


yaw  angle  for  each  location  sampled.  However,  an  interactive  computer  - 
probe  control  system  would  be  needed.  Alternatively,  a method  based  on 
the  approach  developed  in  Appendix  D has  potential  for  success  if  applied 
to  the  proposed  two  probe  system.  The  multiple-sensor  interference 
effects  found  in  the  Dodge  probe,  to  which  the  method  of  Appendix  D was 
applied,  are  not  present  in  the  two  probe  system. 

The  development  of  such  a method  would  require  that  suitable  expres- 
sions be  found  which  would  individually  characterize  the  behavior  of  the 
type  A and  type  B probes  with  respect  to  the  Mach  number,  yaw  angle, 
and  pitch  angle  of  the  flow  to  be  measured.  The  expression  developed  in 
this  study  was  shown  to  characterize  properly  the  behavior  of  cylindrical 
probes  with  respect  to  Mach  number  and  yaw  angle,  but  the  dependence 
on  pitch  angle  was  not  obtained  explicitly. 

Reasonable  agreement  was  rioted  (Table  V-l)  between  the  magnitudes 
of  the  time  averaged  Kulite  probe  impact  pressure  and  the  pneumatically 
averaged  impact  pressure  obtained  using  the  Dodge  probe.  It  was  con- 
cluded that  a Kulite  probe,  with  on-line  calibration  to  an  equivalent  pneu- 
matic probe,  can  be  used  for  quantative  real  time  impact  pressure  measure- 
ments downstream  of  a rotor. 
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CHARACTERISTICS  OF  THE  KULITE-EQUIVALENT  PNEUMATIC  PROBE 
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FIGURE  8.  KULITE  PROBE  CALIBRATION  CURVE 
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FIGURE  9 

EVALUATION  OF  AVERAGING  TECHNIQUES  IN  THE  KULITE  PROBE 

CALIBRATION  TEST 
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FIGURE  10 

PHOTOGRAPH  OF  TYPE  B PNEUMATIC  EQUIVALENT  PROBE 
IN  COMPRESSOR  MOUNTING  APPARATUS 
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FIGURE  1 2d 

BLADE  PAIR  2 KULITE  PRESSURES 
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BLADE  PAIR  2 KULITE  PRESSURES 
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BLADE  PAIR  2 KULITE  PRESSURES 
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FIGURE  I2g 

BLADE  PAIR  2 KULITE  PRESSURES 
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BLADE  PAIR  8 KULITE  PRESSURES 


66 


FIGURE  14.  BLADE  PAIR  2 ZERO  YAW  ANGLES 
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APPENDIX  A 
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DETERMINATION  OF  THE  CHARACTERISTICS  OF 
CYLINDRICAL  IMPACT  PROBES 


A.  1 Calibration  Tests 

A cylindrical  impact  probe  was  constructed  to  have  a tip  with  the 
same  dimensions  as  the  individual  sensors  in  the  Dodge  probe.  The  probe 
is  shown  in  Figure  AI.  Calibration  tests  of  the  probe  were  conducted  in  a 

seven  inch  free  air  jet  at  five  Mach  numbers  varying  from  0.  2 to  0.  5, 

o Q o o 

yaw  angles  from  -90  to  90  and  pitch  angles  from  -40  to  40  . Both 

yaw  and  pitch  angle  were  varied  at  Mach  numbers  0.  2 and  0.  5.  Yaw  angle 
alone  was  varied  at  Mach  numbers  0.  3,  0.  35,  and  0.  4.  A description  of 
the  free  jet  apparatus  is  given  in  Appendix  B of  Ref.  1.  A Prandtl  probe 
was  used  to  measure  the  total  pressure  at  the  same  radial  station  as  the 
test  probe.  A conventional  strain  gauge  transducer  and  scanivalve  were 
used  to  read  the  pressures  from  the  Prandtl  probe,  the  probe  under  cali- 
bration and  the  atmospheric  reference.  Frequent  cross  checks  to  a water 
column  manometer  were  made  to  insure  accuracy.  An  example  of  the 
results  is  shown  in  Fig.  A2, where  the  pressures  measured  in  yaw  surveys 
at  five  Mach  numbers  are  plotted. 

A.  2 Analysis  of  Results 

The  probe  pressure,  Pp,  at  each  yaw  and  pitch  angle  was  reduced 
to  a nondimensional  pressure  coefficient  Cp  defined  as 


ij 
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Cp 


(Al) 


lf/2 Ps  M2 

where  Pt  = reference  impact  pressure  from  the  Prandtl  probe,  Pg  = 
static  pressure  (atmospheric),  M = Mach  number  calculated  using  refer- 
ence impact  pressure  and  atmospheric  static  pressure,  and  = 1,4  = 
ratio  of  specific  heats.  The  values  obtained  for  Cp  were  plotted  against 
the  angle  of  the  flow  to  the  probe  axis.'V  . The  results  are  shown  in 
Fig.  A3.  For  angles  that  were  a combination  of  pitch  and  yaw  angle, 
the  flow  angle  relative  to  the  probe  axis  was  calculated  using  the  geomet- 
rical relationship. 

V = cos"1  (cosO£  cos  0)  (A2) 

where  01  is  the  yaw  angle  and  0 is  the  pitch  angle  of  the  flow  to  the 
cylindrical  axis  of  the  probe.  Table  Al  shows  the  combinations  of  pitch 
and  yaw  angles  from  which  the  angles  in  Figures  A3g  and  A3i  were 
computed. 

The  expression  chosen  to  represent  the  data  at  each  Mach  number 
was  Cp  = A(sin2  B(  - V 0))^.  At  least  squares  algorythm  was  developed 
for  the  computer  which  determined  the  constant  coefficients  K,  B, 
and  N from  given  data.  The  algorithm  is  described  in  Appendix  B. 

Further  examination  revealed  that  for  this  probe  the  Mach  number  depend- 
ence could  be  accounted  for  by  redefining  A = A * M*  so  that  from 
Equation  Al,  the  expression 
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Cp  = A (sin2  B (V-Yo))N  (A3) 

where  Cp  was  defined  as 


and  A,  B,  ^ 0 and  N were  constants,  described  the  behavior  at  all  Mach 

numbers  tested.  This  expression  was  found  to  hold  well  as  long  as  the 

yawed  probe  pressure  did  not  fall  below  the  static  pressure,  which  occurred 

o 

at  a flow  angle  of  approximately  +60  . Application  of  this  expression  (Eq- 
uation A3)  in  a method  for  representing  the  calibration  of  the  Dodge  probe 
is  described  in  Appendix  D.  The  expression  was  used  in  the  present  work 
to  determine  the  yaw  angle  from  measurements  made  with  a cylindrical 
impact  probe  set  at  different  angles  to  the  flow.  This  application  is  dis- 
cussed in  Appendix  B. 

In  the  pitch  angle  surveys  of  the  pneumatic  probe  shown  in  Fig.  A-l 
there  was  a noticeable  assymetry  in  the  distributions  thought  to  be  an  effect 
caused  by  the  probe  stem.  The  assymetry  should  be  noted  in  the  applica- 
tion of  the  third  method  of  calibration  for  multiple  sensor  probes  described 
in  Appendix  D. 
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FIGURE  A-l.  PHOTOGRAPH  OF  0.032"  DIAMETER  IMPACT  PROBE 
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APPENDIX  B 


LEAST  SQUARES  APPROXIMATION  OF  DATA  USING 
Cp  = A(.in2B(  y - Y 0))N 

1.  Derivation  of  A,  B,N  and  Y Q 

From  surveys  of  a cylindrical  impact  probe  in  a known  flow,  pres- 
sure data  is  reduced  to  a set  of  pressure  coefficients  Cpi  and  flow  angles, 
Y For  n sets  of  points  (Cpj,  Y i)  the  coefficients  A,B,N  and  Y0 

form  a least  squares  fit  when  the  expression 
n 

Q=Z  (CprA(sin2B(y.  - H/„))N)2  (Bl) 

i=l 


is  minimized.  A minimum  value  of  Q is  found  when,  simultaneously: 

| [CPi-A(.in2B(  0>>N]  . 

t.in2B(  Yt  - y/ 


<30  = -2  £ 

^A  i=l 


=0  (B2) 


d£_  = 4AN  ^ I (Cpl-A(.ln2B(  Yt  - Y0))N] . [co.BI  Yt  - V „>] 


dQ  = 4ANB 


£ [[Opi-ACin2®,^-^,)^.  [CosB{Yi-Yo)])=0 

[(sin2B(VfrVK0))N’1UsinB(H/i-H/0)]  j 
<HVY0))N].  [(•in2B(Yi-Vo))Nj 

itlV[ln(sin2B(YrY0))] 


(B4) 


dQ  n I [C  .-A(sin2B 
Pl 


(B5) 


0 (B3) 
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The  simultaneous  solution  of  equations  BZ,  B3,  B4  nd  B5  was  achieved 
by  a numerical  procedure  which  sequentially  solved  each  equation  to  find  a 
value  for  the  respective  coefficient  while  the  remaining  coefficients  were 
held  constant.  The  procedure  stepped  through  the  solution  of  each  equation 
in  a cyclic  manner  until  the  simultaneous  solution  for  A,B,  N and 
was  found.  The  solution  of  each  equation  for  the  value  of  the  respective 
coefficient  was  achieved  in  an  iterative  manner  using  Newton's  method  for 
successive  approximations.  The  derivative  used  to  compute  the  iteration 
interval  was  defined  as  the  linear  slope  of  the  function  computed  with  the 
most  recent  pair  of  values  for  the  iterated  coefficient.  The  first  two  values 
of  the  function  were  obtained  using  a fixed  interval  to  increment  the  co- 
efficient. The  solution  was  found  when  the  value  of  the  function  was  less 
than  a specified  small  quantity. 

The  simultaneous  solution  of  the  four  equations  was  obtained  when 
the  value  of  each  respective  partial  derivative  computed  on  the  first 
iteration  step  was  less  than  a specified  small  quantity. 

No  convergence  criteria  was  applied  to  the  sequential  solution 
technique.  After  each  individual  coefficient  was  computed  its  value  was 
held  constant  in  the  calculations  for  the  other  three  coefficients. 

The  BASIC  program  "ITERATIVE  LEAST  SQUARES"  for  this  method 
is  listed  in  Table  B-l.  The  computer  time  required  to  obtain  a simultan- 
eous solution  for  the  coefficients  A,B,N  and  '4/q  depended  on  the  number 
of  data  points  used  and  the  accuracy  of  the  initial  approximation  for  the 
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coefficients.  Experience  gained  in  using  this  procedure  indicated  that  the 
following  guidelines  were  appropriate: 

a.  Critical  data  points  are  in  the  ranges  of  positive  and  nega- 
tive angles  which  define  the  nearly  linear  regions  of  the  curve.  The  three 
nearly  linear  regions  of  the  curve  include,  one  centered  about  =0, 


for  which  is  greater  than  -0.  02,  and  the  two  sides  with  slopes  which 


are  equal  but  of  opposite  sign.  The  two  sides  are  limited  on  one  extreme 
by  the  maximum  and  minimum  angles  and  on  the  other  extreme  by  the 
angles  at  which  significant  curvature  begins.  Thus  six  data  points,  two 
in  each  of  the  nearly  linear  ranges,  are  required  for  a reasonable  fit. 

b.  When  selecting  additional  points  consideration  should  be 
given  to  the  areas  where  the  most  accurate  characterization  of  the  data 
is  required. 

c.  The  initial  estimate  of  should  be  the  angle,  in  degrees, 


at  the  center  of  the  Cp  vs^  curve. 


d.  The  initial  approximation  of  B should  be  1.0. 

e.  The  initial  approximation  of  A should  be  0.  5 less  than  the 


minimum  C of  the  data. 
P 


f.  The  initial  value  of  N should  be  estimated  by  dividing  the 
angular  width,  in  degrees,  of  the  linear  region  for  which  Cp  is  greater 
than  -0.02,  by  18  degrees.  * 


86 


I 


2.  Determination  of  the  Zero  Yaw  Angle 

Using  measurements  of  Cp^  obtained  from  a cylindrical  impact  probe 
in  a known  flow,  the  coefficients  A,B,  N and  are  determined  using 
the  procedure  described  above.  A,B,  and  N now  characterize  the  probe 
with  respect  to  Mach  number  and  the  actual  angle  of  the  flow  to  the  probe 
axis.  If  V is  the  angle  of  the  probe  axis  with  reference  to  some  mech- 
anical scale,  V0  is  a measure  of  the  difference  between  the  zero  on  the 
mechanical  scale  and  the  probe  angle  at  which  the  probe  axis  is  aligned 
with  the  flow  (zero  flow  angle).  This  is  illustrated  in  the  following 
sketch: 


-SB  ~HB  -3B  -20  rlB 


BXI5  or 

FlUH  SWHETRY  — 
(ZERO  FLOW  ANGLE) 


IB  2B  30  HB  SB 


PROBE  ANGLE  SCALE 


When  A,  B,  and  N are  known,  measurements  taken  in  an  unknown 
flow  with  the  same  probe  can  be  used  to  find  which  is  the  angle  of 

the  flow  referenced  to  the  mechanical  scale.  The  pressure  insensitivity 
of  the  cylindrical  probes  near  the  zero  flow  angle  and  the  relationship 
in  Equation  A2  allows  the  assumption  to  be  made  that  for  small  pitch  angles 
the  calculated  zero  flow  angle  of  the  probe  will  be  equal  to  the  zero  yaw 
angle. 
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In  order  to  determine  the  zero  yaw  angles  measurements  of  pressure 
are  made  at  a series  of  probe  angles  about  a reference  value.  Equal 
magnitudes  on  either  side  of  the  reference  are  taken  (i.  e.  = 0, 

♦ Y , + Y rpaV).  The  procedure  to  calculate  the  zero  yaw  angle 

of  the  flow  is  as  follows: 


a.  Reduce  pressure  data  to  values  of  C . In  a known  flow 

P 


this  is  simply  to  use  Equation  A4.  In  an  unknown  flow  the  must  be 


approximated  since  the  total  and  static  pressures  are  not  known.  A suit- 
able approximation  is  obtained  using  the  equation 

N ^ '^max 


Cp=A(sin  B(H/  max)r 


P~;+P. 


P -(fjmnnnun) 


(B6) 


max T 

where;  A,  B and  N are  previously  determined  calibration  constants. 


P.  = ith  impact  pressure  in  the  series 


P = maximum  P.  in  the  series 
max  i 


Pmin  = niinimum  Pj  in  the  series 
Pnrn-n  = next  to  minimum  P.  in  the  series 


'y  = absolute  magnitude  of  the  maximum  and  minimum 
' max 


probe  angles  in  the  survey. 

This  equation  scales  the  pressure  data  to  form  a coefficient  which  has  a 
variation  similar  in  shape  to  the  variation  of  Cp  determined  in  the  calibra- 
tion test. 

b.  Determine  V Q from  the  data  by  solving  equation  B4  for 
V while  keeping  A,  B,  and  N constant.  The  resulting  will  be 
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the  least  squares  solution  for  the  zero  yaw  angle  of  the  flow  with  respect 
to  the  mechanical  probe  angle  scale. 


1 


A BASIC  program  "KAW40"  which  accepts  pressure  data  and  out- 
puts the  value  of  is  listed  in  Table  B2. 
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340  PRINT  "Y0:  "Y0*h 
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APPENDIX  C 

DATA  ACQUISITION  SOFTWARE 

C-l  BASIC  PROGRAM  "DATACQ" 

The  BASIC  program  "DATACQ"  was  written  to  acquire  data  using 
the  Hewlett-Packard  HP  21MX  computer  and  peripherals  described  in 
Section  HI  of  this  paper.  It  is  modular  in  construction  and  can  be  used 
by  a relatively  inexperienced  operator  to  sample  real  time  data  in  a 
variety  of  modes.  When  used  interactively  as  a multipurpose  data  acquisi- 
tion control  program,  a question  and  answer  format  is  used.  For  special- 
ized data  acquisition  tasks,  it  can  accept  by  merging,  user  subprograms 
which  enter  the  modules  (subroutines)  in  "DATACQ"  in  a non-interactive 
manner,  thus  eliminating  the  need  for  keyboard  entries. 

The  key  to  the  flexibility  of  "DATACQ1  is  the  "Selector  Module" 
which  provides  the  control  to  access  all  data  taking  subroutines  on  command. 
A brief  "user  manual"  is  provided  in  this  appendix  to  acquaint  the  user  with 
the  operation  of  "DATACQ".  A program  listing  is  provided  in  Table  C-l. 

C-2  USER  MANUAL  FOR  BASIC  PROGRAM  "DATACQ" 

A limited  amount  of  information  is  provided  in  the  remarks  con- 
tained in  the  program.  This  users  manual  provides  instructions  for  each 
function  preprogrammed  into  "DATACQ".  The  first  step  in  using  "DATACQ" 
is  to  enter  a user  program  that  is  written  separately  then  merged  into  the 
main  program.  The  user  program  must  be  written  between  lines  20  and 
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2000  to  avoid  overwriting  other  program  lines.  The  simplest  user  pro- 
gram that  can  be  used  to  acquire  data  is: 


21  GOSUB  5000 

22  STOP 

This  program  transfers  program  control  to  the  selector  subroutine. 

Subroutine  5000  - Selector  Module 

Upon  calling  subroutine  5000,  the  query 

ENTER  SELECT  VARIABLE 

? 

will  be  displayed.  A number  corresponding  to  the  desired  module  in  the 
following  list  should  be  entered:  * 

1 - Initialize  calibration  coefficient  array 

2 - Enter  calibration  coefficients  from  console 

3 • On  line  probe  calibration 

4 - Synchronized  sampling  of  data 
6 - Free  run  sampling  of  data 

Entry  of  any  number  other  than  the  numbers  listed  here,  0and5,  causes  the 
query  "ENTER  SELECT  VARIABLE"  to  be  repeated.  After  a selected 
module  has  completed  its  operation  "ENTER  SELECT  VARIABLE"  will  be 
repeated  until  a 0 or  a 5 is  entered,  signaling  the  program  execution  to 
return  to  the  user  program. 

Subroutine  4000  - Initialize  Coefficient  Array 

This  subroutine  initializes  the  calibration  coefficients  for  all  data 
channels  to  have  the  values  xQ  = 0 and  X}  = 1.  As  a result,  samples  of 
data  taken  on  all  channels  will  be  displayed  directly  as  voltages. 
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Subroutine  4050  - Enter  Calibration  Coefficients 

This  subroutine  allows  console  entry  of  the  calibration  coefficients 
Xq  and  xj  for  all  16  channels.  Previous  coefficients  for  the  user  selected 
channel  are  lost  upon  entry  of  the  new  coefficients. 

Subroutine  3500  - Time 

If  the  proper  time  of  day  was  entered  into  the  computer  on  startup, 
this  subroutine  returns  with  Rl,  T2,  and  T3  set  to  the  time  of  day  in  hours, 
minutes  past  the  hour,  and  seconds  past  the  minute  respectively.  T is 
set  to  the  time  of  day  in  seconds. 

Subroutine  3600  - Time  Scheduling 

Caution  - This  subroutine  may  be  called  only  once  during  "DATACQ" 
execution.  A second  call  will  result  in  an  indefinite  waiting  state  in  pro- 
gram execution  requiring  an  "abort"  command  to  return  computer  control 
to  the  operator. 

The  purpose  of  this  subroutine  is  to  allow  repetition  of  a task  with 
start  time,  number  of  repetitions,  and  the  interval  between  repetitions 
defined  by  the  user.  The  task  to  be  repeated  must  be  programmed  with 
the  first  program  line  number  being  100,  300,  500,  or  900  and  a "Return" 
command  as  the  last  program  line.  Before  calling  the  subroutine,  the 
variable  V9  must  be  set  to  a value  equal  to  the  first  line  number  of  the 
task  to  be  repeated. 

On  completion  of  the  required  number  of  task  repetitions,  execution 
of  a users  program  will  continue  with  program  line  30. 
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Caution  - The  user  must  insure  that  the  time  interval  entered  is 
sufficient  to  complete  the  task  scheduled.  Overlapping  of  scheduled  tasks 
is  not  allowed. 

Subroutine  2000  - RPACE  Subroutine 

This  subroutine  is  used  for  synchronized  sampling  of  a selected  data 
channel.  Before  calling  the  subroutine  the  following  variables  must  be 
defined. 

A1  = 32768  + 256  * (Blade  Pair)  + SAMPLE  LOCATION  (0-255) 

N1  = number  of  data  samples  to  be  taken  at  the  specified  location 

(Maximum  value  2000) 

N2  = number  of  channel  to  be  sampled 

The  returned  variables  and  their  values  are: 

A2  = (15x106)/RPM 

C(20,  100)  = Data.  The  first  N1  elements  of  this  array  (in  row 
order)  contain  the  values  of  the  data  samples. 

C0  = The  numerical  average  of  N1  samples 

Subroutine  2500  - Free  run  sampling 

This  subroutine  is  used  for  non- synchronized,  "Free  Run"  sampling 
of  a selected  channel.  Input  variables  N1  and  N2  are  defined  as  in  the 

RPACE  subroutine.  The  N1  samples  requested  are  taken  at  10  micro- 

♦ 

second  intervals.  Output  variables  C(20,  100)  and  C0  are  as  defined  in 


the  RPACE  subroutine. 
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Subroutine  7300  - Subroutine  to  sample  pressure 

This  routine  will  sample  channels  designated  as  pressure  channels 
and  a channel  designated  as  a reference  pressure  channel  in  the  free  run 
mode,and  reduce  the  voltage  data  on-line  to  pressures  using  the  expres- 
sions 

Pr  = <x0)r  + (*1)r*(E)r 

and 

Pp  = (x0)p  4 (Xl)p‘  (E)p  + Pr 

where  the  subscript,  r,  refers  to  the  reference  channel  and  subscript  p, 
refers  to  the  pressure  channel.  xQ  and  Xj  are  the  calibration  coefficients 
for  the  selected  channel,  E is  the  ensemble  average  of  the  sampled  volt- 
ages, Pr  is  the  reference  pressure  and  Pp  is  the  calculated  pressure  for 
the  selected  channel. 

Inputs  are  defined  interactively  by  answering  the  queries  printed 
on  the  CRT  console  or,  for  non-interactive  use,  values  may  be  pre- 
assigned to  the  listed  variables  and  the  subroutine  entered  at  line  7340. 
The  variables  which  must  be  defined  in  the  non -inter active  mode  are: 

E3 -Channel  # for  the  reference  pressure  or  enter  16  for  input  at 
the  console 

E4-Number  of  free  run  samples  desired  on  reference  pressure 
channel 

E5-Channel  number  of  the  pressure  channel  to  be  sampled 

E6-Number  of  free  run  samples  desired  of  the  pressure  channel 
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The  subroutine  prints  the  pressure,  reference  pressure,  and  pres- 
sure channel  average  voltage  on  the  console  CRT. 

C-3  CONTROL  PROGRAM  FOR  COMPRESSOR  TEST 

A program  was  written  to  be  merged  into.  "DATACQ"  which  caused 
data  to  be  output  on  the  tape  punch  and  sent  to  the  HP  9830  calculator.  The 
program  was  used  for  the  compressor  test  described  in  Section  V.  A list- 
ing of  the  control  program  is  given  in  Table  C-2. 
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3500  REH  TIME  SAMPLE  SUBROUTINE  RETURNS  WITH 
3)10  REM  TUMOURS, T2:HINUTES,T3=SEC0NDS 
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7500  REM  SUBROUTINE  FOR  PACED  SAMPLING 

7501  REM  AVG  VOLTAGE  FOR  EACH  DISCRETE  SAMPLE  POINT  IN  ARRAY  D(POINT) 

7502  REM  INTERACTIVE  ENTRY  SUB  7500 

7503  REM  PREDEFINED  ENTRY  7540|J2=dL  PR,J3=SAM  SP,J4=1  PER  REV.N1.N2 
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APPENDIX  D 


DERIVATION  OF  MULTIPLE  SENSOR  PROBE  CALIBRATION 
USING  SINGLE  IMPACT  PROBE  CHARACTERISTICS 


Using  the  characteristics  of  the  cylindrical  impact  probe  discussed 
in  Appendix  A,  which  was  similar  in  geometry  to  the  individual  sensors  of 
the  Dodge  probe,  a calibration  was  derived  for  the  Dodge  probe. 

The  following  definitions  were  used: 

P^  - Pressure  at  the  Dodge  Pj  sensor 

P - Average  of  the  Dodge  P,  and  P,  sensor  pressures 
23  4 

P^  - Pressure  at  the  Dodge  P^  sensor 

V|/.  - The  angle  of  the  flow  with  respect  to  the  sensor  axis 

0£  - Yaw  angle  with  respect  to  axial  direction 

0 - Pitch  angle  with  respect  to  the  plane  normal  to  the  probe 

shaft 

M - Mach  number 

Pg  - Static  pressure 

P - Total  pressure 

c 

y - Ratio  of  specific  heats  ( = 1.4  was  used) 

In  Appendix  A it  was  shown  that,  for  a single  impact  tube,  the 
pressure  coefficient  defined  as 


C = 
P 


Pi-Pt 

i p.™‘ 


was  closely  approximated  by  the  expression 
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r 


i 


cp  = A [«in2B(  yr  yo)]N  (2) 

where  the  flow  angle  ^ was  determined  from  pitch  and  yaw  angle  by 
M^cos’VcosQCiCOS  0j)  (3) 

For  the  Dodge  probe,  when  rotated  to  balance  the  pressures  at  sensors 
P^  and  P^,  the  following  table  lists  the  values  of  the  angles  for  each  of 
the  sensors 


Dodge 

sensor 

1 

2-3 

4 

«, 

0 

(X 

0 

fli 

0 

0 

0-* 

% 

0 cos 

‘‘(CO.®  COS0) 

0-0 

where 

OC  = angle  between  axis  of  sensor  2 or  3 and  sensor  1 

7 = angle  between  axes  of  sensor  4 and  sensor  1 


The  quantities  which  have  been  used  previously  to  represent  the  probe 

I 

characteristics  are: 


P 


1 


P 


1 


and 


0- 


PUPJ3 

“PT* 


Using  equation  (2)  in  equation  (4) 


(4) 


(5) 


cp.  - °p« 


[•in2B  VP  - [iln2B  V,]N  • 
[.U2B  V,1N-[.i«2B  V23]N 
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Using  equation  (3), 


[sin2Ba]N  - [sin2B(0-0)]N  ^ 

[sin2B0]^-[*in2B(cos" *(cos  (X  cos  0)>  )]^ 

Note  that  if  the  pressure  coefficient,  Cp  is  independent  of  Mach  number, 
then  equation  (6)  gives  pitch  angle  (0)  in  terms  of  sensor  measurements  ( ). 

There  must  be  no  interference  between  sensors,  however  for  equation  6 
to  describe  the  probe  characteristics  correctly.  Using  Equation  (5) 


Q = — P1  ~ CP23  s A[sin2Bg]N  - A [sin2B(cos“1  ( cos  (X  cos  p) 

7 Si+^ r A[*ia2Bg]N  - [l+ii-^r 

1P.M2  /M  2 


A[ain2B0r  - -£-;z  [l  + %’J  M2]" 


PSM‘ 


Defining 


X = 2 


7m2 


[ 1 + Izl  M2]^  1 


equation  (7)  may  be  written 

X = — “1)  (sin2  (B0)  ) - JL  (*in2B(cos  * (coaOC  cos  t)  ))  ] 

P JJ 

If  the  pitch  angle  is  known,  equation  (9)  gives  the  Mach  number,  M,  for  a 
measured  value  of  f2. 

Alternatively  if  we  define 


6 -ft  7 - 


Cp  1 ■ CP< 


cpi  + x 


then  equation  (9)  becomes: 


= A -1)  (sin2  (BO)  ) - -g— (sin2B  ( 0-0)  )N]  (II] 


Equations  (6),  (8),  (9)  and  (11)  were  used  to  derive  a calibration 
for  the  Dodge  probe.  The  coefficients  A,  B,  and  N which  represent  the 
characteristics  of  the  single  impact  tube  similar  to  the  sensors  of  the 
Dodge  probe  were  determined  as  in  Appendix  B. 

From  measurements  taken  in  a known  flow  at  zero  pitch  angle 
(0  = 0)  with  the  Dodge  probe  the  quantities  J2 , If  and  (5  were  calculated 
using  equations  4,  5 and  10. 

The  angle  (X  was  computed  from  equation  9. 

First,  x was  calculated  from  equation  (8)  using  the  measured 
Mach  numbers.  Then  from  equation  (9), 

_ . At.in2  BCf  )N 

P 

or  , _! 

a r 


The  angle  7 was  computed  similarly  from  equation  11: 

if/  X ^ , vJ-  1 1 


? = sin 


■)  2N  ] 


Having  established  Q and  J from  measurements  the  procedure  to  com- 
pute velocity  using  values  of  and  & measured  in  an  un- 


known flow  was: 


1)  Assume  0=0 


2)  Calculate  the  right  hand  side  of  equation  6 = ^ 

3)  If  \f-7]  ~ 0 then  a solution  for  0 was  found.  Proceed 


to  calculate  M. 
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4) 


Calculate 


&7/&0 


5)  Set  0 = 0 4 

6)  Return  to  step  2 with  the  new  value  of  0 

Repeat  the  iteration  .until  a solution  for  0 was  found 
To  calculate  M: 

7)  Calculate  the  right  hand  sides  of  the  equations  (9)  and  (11). 

Compute  the  average  = x . 

8)  Assume  a starting  value  for  M 

9)  Calculate  x using  equation  (8) 

10)  If  |x  “ 3(j  0 then  a solution  was  found  and  the 

unknown  flow  was  determined 

11)  Calculate  d */  dM 

X - 7 

12)  Set  M = M + dydM 

13)  Return  to  step  9 with  the  new  value  of  M 

Repeat  the  iteration  until  a solution  for  M was  found. 

This  method  (method  III)  for  calibration  of  the  Dodge  probe  was  tested  using 
data  previously  obtained  with  the  Dodge  probe  in  a known  flow.  Calibration 
methods  I and  II  were  used  to  reduce  the  same  data  to  pitch  angle  and 
velocity  magnitude.  (Each  data  point  corresponded  to  values  of  pitch 
angle  and  velocity  magnitude  which  were  controlled  and  known  during  the 
calibration  test.)  Table  D-l  is  a comparison  of  the  accuracies  achieved 
using  the  three  methods. 
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It  can  be  seen  that  method  III  was  the  least  accurate  method  of 
representing  the  calibration  of  the  Dodge  probe.  However,  the  mutual 
interference  of  the  probe  tips  is  probably  the  cause  of  the  large  error.  In 
the  2 probe  system  of  velocity  measurement  there  can  be  no  interference 
between  probes  and  hence  this  method  should  be  further  considered  for 
that  application. 
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MACH  NO.  PITCH  ANGLE 


CALIBRATION 

AVG. 

MAX 

AVG. 

MAX 

MACH  # 

METHOD 

ERR  % 

ERR  % 

ERR  DEG 

ERR  DEG 

.35 

I 

2.6 

3.2 

0.5 

-0.9 

.35 

II 

0.5 

-1.2 

0.6 

-1.0 

.35 

III 

4.  5 

7.3 

0.8 

-1.  1 

.44 

I 

1.  5 

-2.2 

0.4 

-1.0 

.44 

II 

0.3 

-0.  6 

0.4 

-0.  6 

.44 

in 

4.7 

8.3 

1.0 

-1.4 

.49 

I 

2.0 

-4.2 

0.6 

-1.8 

.49 

n 

0.  5 

0.  9 

0.4 

0.  7 

.49 

iii 

4.9 

8.4 

1.4 

-1.  6 

Average  and  maximum  values  of  the  errors  are  shown  for  the  five  pitch 
angle  settings  -5°,  0°,  5°,  10°  <&  15  . 


TABLE  D-l  SUMMARY  OF  ERRORS  IN  THE  USE  OF 
THREE  CALIBRATION  METHODS  FOR  THE  DODGE  PROBE 
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